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ABSTRACT 

We perform Monte Carlo simulations of cosmic ray-induced hard X-ray radiation from 
the Earth's atmosphere. We find that the shape of the spectrum emergent from the 
atmosphere in the energy range 25-300 keV is mainly determined by Compton scat- 
terings and photoabsorption, and is almost insensitive to the incident cosmic-ray spec- 
trum. We provide a fitting formula for the hard X-ray surface brightness of the atmo- 
sphere as would be measured by a satellite-born instrument, as a function of energy, 
solar modulation level, geomagnetic cutoff rigidity and zenith angle. A recent mea- 
surement by the INTEGRAL observatory of the atmospheric hard X-ray flux during 
the occultation of the cosmic X-ray background by the Earth agrees with our pre- 
diction within 10%. This suggests that Earth observations could be used for in-orbit 
calibration of future hard X-ray telescopes. We also demonstrate that the hard X-ray 
spectra generated by cosmic rays in the crusts of the Moon, Mars and Mercury should 
be significantly different from that emitted by the Earth's atmosphere. 
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Starting in the 1960's, observations from bal- 



loons (iPeterson et al.l Il973l : iDaniel fc Stephens! 1 1974 



19741: ISchwartz fc PetersonI 



Shaw et al.l 



Willett fc Mahonev 
20031 ) 



19741: llmhof et all 


197e; 


Share & Murohvl 


2001; 



have shown that the Earth's atmo- 
sphere is bright in hard X-rays and gamma-rays owing to 
its bombardment by cosmic rays of Galactic origin. This 
is the case despite the fact that the energy deposited by 
the cosmic rays into the atmosphere is some nine orders of 
magnitude smaller than the total energy budget coming to 
us from the Sun and that only ~ 1% of the energy deposited 
by the cosmic rays is re-radiated into space as gamma-rays 
and much less as hard X-rays. 

Atmospheric hard X-ray emission is of interest not only 
because it provides insight into the interaction of energetic 
cosmic particles with air, but also in connection with the 
possibility of measuring the spectrum of the cosmic X-ray 
background (CXB) using the Earth as a shield. In January- 
February 2006 a series of observations following this ap- 
proach were performed by th e INTEGRAL orbital obser- 
vatory l|Churazov et al.l |2007| ) and drew our attention to 
the problem of atmospheric X-ray radiation. In observa- 
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tions of this kind, cosmic ray-induced emission on the one 
hand causes significant contamination of the CXB decre- 
ment (due to occultation by the Earth) at energies ~ 30- 
100 keV, but on the other, could in principle be used for 
calibrating a detector's response in the high-energy region 
(substantially above 100 keV) where atmospheric emission 
dominates over the CXB signal. During the INTEGRAL ob- 
servations the Earth covered ~ 100 sq. deg of the sky and 
was as b right as ~ 200 mCra b in the spectral band 100- 
200 keV jChurazov et al.ll2007l ). 

Despite the long hist ory of experimental (see references 
abov e ) and theo r etical (IPus kin '1970*; 'Daniel fc Stephend 
1 19741 : iLind Il975l : ICraser fc S chonfcldcr 1977; Dean ct 4} 
'1989^) studying, there remains significant uncertainty re- 
garding the spectrum of atmospheric X-ray and gamma-ray 
radiation. However, the nuclear and electromagnetic pro- 
cesses responsible for the generation and propagation of 
X- and gamma-rays in the atmosphere, although complex, 
are well understood and already incorporated in computer 
codes. Also the incident energy spectrum of Galactic cos- 
mic rays has been measured with high precision (see e.g. 
iGaisser et al.ll200ll for a review) and there is a good under- 
standing of how it is modified below ~ 10 GeV by the mag- 
netic fields of the Solar system and Earth. All this suggests 
that the energy and angular distributions of atmospheric X- 
ray and gamma-ray emission should be predictable by sim- 
ulations with fairly good (~ 10%) accuracy. 
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We have performed suc h a numerical modellin K using 
the popular toolkit GEANT llAgostinelU et al.ll20oj ). focus- 
ing on making predictions for observations of the Earth from 
space in the hard X-ray band 25-300 keV. We show below 
that the hard X-ray emission of the Earth is characterized 
by a well-defined spectral shape (Sect.[2| and that its inten- 
sity depends in a predictable way on the phase of the solar 
cycle and on the geomagnetic rigidity map corresponding 
to a given observation (Sect.[31). Therefore the Earth could 
become a good absolute calibrator for future orbital hard X- 
ray telescopes, in particular when used in combination with 
accurate monitoring of the cosmic-ray spectrum by large- 
area, high-resolution particle physics detectors such as the 
Alpha Magnetic Spe ctrometer to oper ate on the Interna- 
tional Space Station (|Barao et al.ll2003 V 

In addition we demonstrate (Sect. |3| that the cosmic 
ray-irradiated solid crusts of the Moon, Mars and Mercury 
should generate hard X-ray spectra different from the one 
produced in the Earth's atmosphere. 



2 SIMULATIONS 

Our Monte Carlo code is based on the GEANT4 software 
package (release 8.1), which was developed at CERN as an 
extension of GEANT3 and ha s many applications a cross 
particle physics and astronomy (|Agostinelli et al"]|2003l ). We 
adopted the physics model originally constructed for eval- 
uating the effects of exposing the test masses of the Laser 
Interferometer Spac e Antenna to energetic particles in space 
l|Arauio et al.ll2005l ). This model includes a state-of-the-art 
description (in particular, intra-nuclear cascade models) of 
hadronic interactions of cosmic-ray protons and a-particles, 
and their air-shower products over many orders of magni- 
tude in energy, as well as the low-energy electromagnetic 
package in GEANT4, which extends the standard electro- 
magnetic package by carefully treating such processes as 
Rayleigh scattering, photoabsorption and fluorescence, im- 
portant in the X-ray regime. 

Unless otherwise stated, we use a plane-parallel ap- 
proximation for the Earth's atmosphere. Specifically, our 
adopted geometry is a slab of height 96 km and of very 
large (2 x 10'' km) width and length. The air consists of 
atoms of nitrogen, oxygen and argon with mass fractions of 
75.5%, 23.1% and 1.3%, respectively. The air's density de- 
creases exponentially with altitude from a sea-level value of 
1.2 mg cm~'^ with a scale height of 8 km. The total column 
density of our model atmosphere is thus ~ 960 g cm^^. 
This model is designed to approximate the real structure 
of the Earth's atmospher^B below ~ 100 km. The remain- 
ing Earth's atmosphere above this level is so rarefied (~ 
1 mg cm~^) that only a tiny fraction (< 10~^) of the total X- 
ray and gamma-ray radiation escaping from the atmosphere 
into space is produced there, since the bulk of the emergent 
flux is produced in the upper several tens of g cm~^ of the 
atmosphere (above ~ 20 km from the sea level). 

Importantly, by performing many test runs, we found 
that the properties of the emergent hard X-ray radiation 
remain unchanged within the statistical uncertainties if we 

^ See e.g. |http: / /modeIweb.gsfc.nasa.gov| 



increase or decrease severalfold the total height, scale height, 
or the sea-level density of the atmosphere, as long as the to- 
tal column density of the atmosphere remains higher than 
~ 200 g cm"'^. It is the chemical composition of the air that 
largely determines the outcome of simulations. Furthemore, 
we also carried out a series of simulations using a realistic 
spherical geometry for the atmosphere and convinced out- 
selves that the plane-parallel approximation is adequate for 
treating the problem under consideration (see Sect. 12. Il be- 
low). 

We emphasize that we use our simulations to formu- 
late predictions specifically for observations of the Earth 
from satellite orbits. This implies that the presented results 
should not be applied to observations performed from bal- 
loon (i.e. intra-atmospheric) altitudes. 

2.1 Spectra induced by monoenergetic protons 

We first performed a series of simulations in which the at- 
mosphere was bombarded by monoenergetic protons. The 
protons were assumed to come from a uniform background 
subtending the 2n space above the atmosphere. This im- 
plies that the flux density of incoming particles dF/dfi <x /x, 
where fi — cosO and 9 is the zenith angle. We accumulated 
the energy spectra (above an imposed threshold of 1 keV) 
and angular distributions of X- and gamma-rays produced in 
the proton-induced cascades and emergent from the atmo- 
sphere. Similar information was collected for the emergent 
electrons and positrons, to later assess the role of re-entrant 
albedo particles. 

We typically sampled (50,000-125,000)/£'p protons in 
a single run, where Ep (GeV) is the proton kinetic energy. 
By performing a number of test runs, we convinced our- 
selves that this number of sampled particles enables recov- 
ering outgoing hard X-ray to gamma-ray spectra (between 
~ 25 keV and ~ 10 MeV) with a bin-to-bin statistical noise 
of less than 5% when a logarithmic grid of 10 energy bins per 
dex is used, and with a scatter in the amplitudes of spec- 
tra obtained in different runs of less than 3%. These lev- 
els of statistical uncertainty characterize all spectral plots 
presented below. We note that the accuracy of our simula- 
tions becomes signiflcantly worse at photon energies below 
~ 25 keV and above ~ 10 MeV (see in particular Fig. [3] 
below), in the regions lying outside the focus of this paper. 

Figures [T] and [5] show the emergent zenith angle- 
integrated photon spectra for Ep = 1 and 10 GeV. One 
can see that several photon production mechanisms provide 
comparable contributions to the outgoing spectra. These in- 
clude: i) free-free emission by electrons and positrons (dis- 
tributed over several orders of magnitude in photon en- 
ergy), ii) electron-positron annihilation (mostly in a line at 
511 keV) and iii) high-energy processes including nuclear de- 
excitation (in the ^ 1-15 MeV rang^and tt" decay (a broad 
bump near 100 MeV). The energy albedo of the atmosphere 
(ratio of the emergent radiation flux over the incident energy 
flux of cosmic rays) ranges from ~ 1% for 1 GeV protons to 
- 0.1% for 250 GeV protons. 

^ Note that because of the short lifetimes (typically a tiny frac- 
tion of a second) of the excited state, atmospheric gamma -ray 
lines may be considered prompt (see e.g. iRamatv et al]ll979l .l 
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Figure 1. Top: Spectral energy distribution (in E'^dN-y/dE 
units) of photons emergent from the atmosphere exposed from the 
top to an isotropic distribution of protons of energy Ep = 1 GeV. 
The thick solid line shows the total spectrum while the other lines 
show its different components according to the photon production 
process: free-free emission (solid), electron-positron annihilation 
(dashed), other processes including nuclear deexcitation and tt" 
decay (dash-dotted). Also shown are the spectra (E^dN^/dEc) of 
secondary electrons and positrons escaping from the atmosphere 
(dotted). Bottom: The same but for Ep = 10 GeV. 



As is evident from Fig. [5J most of the escaping photons 
with final energies below ~ 511 keV have experienced sin- 
gle or multiple Compton down-scatterings and thus lost a 
significant or dominant fraction of their initial energy due 
to electron recoils. The emergent hard X-ray spectrum has 
thus lost much of the memory about its parent emission pro- 
cesses and as a result is fairly similar to the Comptonized 
spectra appearing in various astrophysical situations where 
hard X-rays or gamma-rays propagat e through a dense and 
cold environment , e.g. in supernovae (|Grebenev fc SunvaevI 
1 19871 ). I n particular, su ch a spectrum was observed from SN 
1987A (jSunvaevet aLlflgST.! . 



2.1.1 Adequacy of the adopted atmosphere model 

Our using of a plane-parallel model for the Earth's atmo- 
sphere is motivated by the fact that all of typical free paths 
of particles and photons in the problem under consideration 
are less than a few ten km and are thus much shorter than 



Figure 2. Top: Spectral energy distribution (in E'^dN.y /dE units, 
thick solid line) of photons emergent from the atmosphere exposed 
from the top to an isotropic distribution of protons of energy 
Ep = 1 GeV. The thin line shows the spectrum of the initial 
energies (i.e. before Comptonization) of the escaping photons. 
Bottom: The same but for Ep = 10 GeV. 



the Earth radius. Nevertheless, to verify the adequacy of 
the adopted approximation we performed a series of more 
realistic simulations, in which the Earth's atmosphere was 
modeled as a spherical shell with an inner radius of 6400 km 
and an outer radius of 6496 km. The vertical structure and 
chemical composition of the atmosphere were adopted ex- 
actly the same as in our reference (slab-geometry) model. 
In Fig. Owe compare the angular-integrated emergent spec- 
trum produced by 1 GeV protons in the spherical atmo- 
sphere with that produced in the plane-parallel one. As ex- 
pected, there is no significant difference between these two 
spectra apart from the statistical noise. We found this to be 
the case also for other incident proton energies. These results 
confirm that the plane-parallel approximation is appropriate 
for our study (see also Sect. 12. 0] below). 

We also show in Fig. [3] a simulated spectrum for a 
plane-parallel atmosphere that has the same sea-level den- 
sity (1.2 mg cm~3) and scale-height (8 km) as our reference 
model but is only 48 km (instead of 96 km) high. In this 
case too there are no significant deviations in the outgoing 
spectrum relative to the one computed within our standard 
model. We present this example to demonstrate the very 
weak sensitivity of the emergent spectrum to changes in the 
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Figure 3. Top; Spectral energy distribution (in E^dN-y/dE 
units) of photons emergent from the atmosphere exposed from the 
top to an isotropic distribution of protons of energy Ep = 1 GcV 
for different models of the atmosphere; a slab geometry with a 
total height of 96 km (the reference model) - the dashed line, 
a spherical shell geometry with an inner radius of 6400 km and 
the same vertical structure as the reference model - the dashed 
line, and a slab geometry with a total height of 48 km and the 
same vertical structure as the lower half of the reference model - 
the dash-dotted line. Bottom; The ratio of the spectra obtained 
for the spherical and 48-km slab models to that obatined for the 
reference model — the solid and dash-dotted line, respectively. 



setup of the atmosphere other than its chemical composi- 
tion, which was mentioned at the beginning of Sect. [2] 

In addition, Fig. [3] provides a sense of the statistical 
uncertainties pertaining to our simulations, which, as we 
have already mentioned, are at the level of a few per cent 
over the 25 keV-10 MeV energy range. 



2.1.2 Fitting formulae for the emergent hard X-ray 
spectrum 

On performing many simulations, we proved (see Fig. 2]) 
that the emergent spectrum between 25 and 300 keV is 
barely sensitive to the energy of the parent cosmic proton 
in the range of interest to us, namely from ~ 400 MeV to 
~ 250 GeV (see Sect. |3}. The outgoing photon spectrum 
over 25-300 keV can be well approximated by the following 
fitting formula: 



C 



dNj_ 

1e ^ (S/44 keV)-5 + [E/U keVy - 



keV~ 



(1) 



where the coefficient C depends on the primary particle en- 
ergy as is discussed below. 

The photon spectrum (Fig.|4]| peaks at ~ 50-60 keV. At 
energies below 40 keV the spectrum exhibits a rapid decline 



Figure 4. Emergent spectra (solid lines, note the different units 
compared to the previous figures) in the X-ray to soft gamma-ray 
range for different energies of incoming cosmic protons (labels 
next to the curves). The dashed lines show the analytic approx- 
imation for the 25—300 keV range given by equation ((TJ. This 
approximation breaks down at Ep <400 MeV. 



(oc E^) due to Compton down-scatterings. Above ~ 300- 
400 keV, where the escaping photons have experienced at 
most few Compton scatterings, the spectrum (including the 
strength of the 511 keV line) starts to exhibit a noticeable 
dependence on the incident proton energy. This is also the 
case for photons with energies less than ~ 25 keV, most of 
which were born via bremsstrahlung in the skin of the at- 
mosphere with surface density ^ several g cm~^ and whose 
spectrum retains signatures of the energy distribution of the 
parent secondary electrons. We note that at these low ener- 
gies the emergent X-ray flux is very weak due to photoab- 
sorption. 

We next address the amplitude of the resulting photon 
spectra. In Fig. [S]the coefflcient C appearing in equation ((T} 
is shown as a function of the incident proton energy Ep . The 
results of simulations in the range 400 MeV-250 GeV are 
well described by the following expression: 

r (F) = (9) 

(Sp/0.6 GeV)-i-8 + (Ep/O.G GeV)-o-7 ' ^ ' 

which allows one to calculate the emergent intensity for a 
given incident spectrum of cosmic protons, as we do in fjS] 
below. 



2.1.3 Angular dependence of emergent radiation 

We proceed to compare spectra emergent at different zenith 
angles. We found (see Fig. |6]) that the spectral shape be- 
tween 25 and 300 keV remains stable and as given by equa- 
tion IT]), while the flux (Fig. [T]) approximately obeys the law 
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Figure 5. Normalization of the emergent 25—300 keV spectrum, 
equation ([TJ, as a function of proton energy Ep. Shown are values 
derived from simulations (points) and the fitting curve (solid) 
given by equation ([2j. An arbitrarily scaled linear dependence is 
shown for comparison (dotted). 
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Figure 6. Top: Spectra for Ep = 1 GeV, emergent at dilfer- 
ent zenith angles {fi = cos 6): fi =0-0.2 (solid), 0.2-0.4 (dotted), 
0.4-0.6 (short-dashed), 0.6-0.8 (long-dashed) and 0.8-1.0 (dash- 
dotted). The spectra were divided by /i(l-|-/i). Bottom: The same 
but for Ep = 10 GeV. 



dF oc /i(l + ^)d/i. 



(3) 



This means that the outgoing hard X-ray radiation ex- 
hibits some central brightening as compared to the classical 
Lambert case, dF oc /id^, corresponding to black-body radi- 
ation. The angular distribution of the atmospheric hard X- 
ray emission is also different from the law of darkening per- 
taining to the other classical case of a semi-infinite scatter- 
ing atmosphere wi th radiation coming from deep inte riors: 
dF cx ^.{l+2n)dn l|Chandrasekha j| 1960l : ISobolevl(l963h . The 
resulting distribution actually resembles that of the Comp- 
tonized radiation escaping from a purely scattering disk 
with optical thickness tt ~ a few (Titarchu k fc SunvaevI 
This refiects the fact that although cosmic rays gen- 
erate radiation over several hundred g cm~^ deep into the 
atmosphere, only radiation produced in the top few ten 
g cm~^ of air can efficiently survive multiple Compton down- 
scatterings (as well as photoabsorption) and escape into 
space in the form of hard X-rays. 

On the contrary, the emergent gamma-rays (with ener- 
gies ^ 1-2 MeV) exhibit some limb brightening (enhanced 
flux at small /i's and reduced flux at large /i's), see the bot- 
tom panel of Fig. [7]), which can be attributed to the fact 
that escaping gamma quanta have experienced at most a 
few Compton scatterings. 

In Fig.[7]we also compare the result obtained using our 
standard atmosphere model with one computed using the 
spherical setup of the atmosphere, described in Sect. [2TTTT] 
Similarly to the angular-integrated spectrum, there are no 
significant differences between the two derived angular de- 
pendences of the emergent fiux. Specifically, the difference 
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Figure 7. Emergent flux for Ep = 10 GeV as a function of the 
zenith angle {fi = cos 6) in different photon energy ranges: 25— 
100 keV, 100-500 keV, 500 keV-2 MeV and above 2 MeV. The 
solid and dashed lines show the results obtained using the plane- 
parallel and spherical atmosphere models. In the upper two panels 
also a well-matching ~ ^(1 -|- fj.) dependence is shown (dotted). 
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Figure 8. Emergent spectra in tlie X-ray to soft gamma-ray range 
generated by o-particles of energy Ea = 4 and 32 GeV (solid 
lines) compared to the spectra generated by protons of energy 
Ep = Ea (dashed lines) and by four protons each of energy Ep = 
Ea/4: (dotted lines). 



between the two approximations is less than 3% of the aver- 
age flux, (dF/dfi), in each bin including the < ^ < 0.2 one 
(corresponding to zenith angles between 78.5° and 90°). We 
do expect though the plane-parallel approximation to even- 
tually break down at grazing angles approaching 90°. How- 
ever, at these angles only a tiny fraction of the total flux 
from the Earth is emitted (for example the fraction of the 
total flux emitted at /i < 0.2 is ~ 3%), which implies that 
our plane-parallel model should be applicable to virtually 
any space-born experiment. 



Figure 9. Emergent spectra (solid lines) in the X-ray to soft 
gamma-ray range generated by cosmic electrons of different en- 
ergies. The dashed lines show the analytic approximation in the 
25—300 keV range given by equation JTJ. 



2.3 Cosmic electrons and positrons 

A significant fraction of Galactic cosmic rays consists of 
electrons (with a small fraction of positrons) . We therefore 
also performed simulations for monoenergetic electrons. The 
resulting hard X-ray spectra are not significantly different 
from the case of protons and in the energy band 25-300 keV 
equation remains a good approximation of the spectra 
(Fig. m. We also tested positrons and, as expected, obtained 
nearly identical results. 

In Fig. [To] the spectral normalization Ce - to be substi- 
tuted for C in equation (IT| - is shown as a function of the 
electron energy Ec. This dependence can be approximated 



2.2 Alpha particles 

A major component of Galactic cosmic rays is He nu- 
clei. GEANT4 permits to handle inelastic scatterings of a- 
particles up to Ea ~ 40 GeV. We therefore ran a number of 
simulations of bombardment of the atmosphere by monoen- 
ergetic Q-particles and found (see Fig. [Si that the resulting 
hard X-ray spectra are nearly identical to those induced by 
protons of the same energy as the Q-particles {Ep = Ea)- If 
the spectrum generated by an a-particle of energy Ea is in- 
stead compared with that generated by four protons each of 
energy Ep — Ea/'l (a possibility suggested by the simplest 
superposition theory), then the match becomes somewhat 
worse. We can therefore extend the use of equation ((5)), orig- 
inally derived for protons, to a-particles, in which case one 
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(4) 



Similarly to the case of protons, in iJSlbelow we integrate 
equation (|4]) over the energy spectrum of cosmic electrons 
bombarding the Earth. 



2.4 Re-entrant albedo particles 

So far we have implicitly assumed that only energetic par- 
ticles coming from space can induce X-ray and gamma-ray 
emission in the atmo sphere. However, already early stud- 
ies (e.g. IPuskinI IT970l ) indicated that additional radiation 
should be generated by re-entrant albedo particles, i.e. sec- 
ondary electrons and positrons that were born in air show- 
ers, escaped from the atmosphere but were then deflected 
back into it by the geomagnetic fleld. Although an accurate 
evaluation of this effect would require more detailed sim- 
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Figure 10. Normalization of the emergent 25—300 keV spectrum, 
equation |(TJ , as a function of electron energy Ee ■ Shown are values 
derived from simulations (points) and the fitting curve (solid) 
given by equation ((4j. An arbitrarily scaled linear dependence is 
shown for comparison (dotted). 



Figure 11. Top; Emergent spectra produced by the re-entrant 
albedo electrons (dotted) and positrons (dashed) generated by 
protons of energy Ep = 1 GeV. The spectrum produced during 
the primary cascade, multiplied by 0.04, is shown for comparison 
(solid line). Bottom: The same but for Ep = 10 GeV. 



ulations involving magnetic fields, it is possible to make a 
first-order estimate. For this purpose, we used the recorded 
spectra of secondary electrons and positrons induced by cos- 
mic protons and emergent from the atmosphere (Fig. [1]) as 
new input spectra for bombarding the atmosphere. This is a 
reasonable procedure since the low magnetic rigidities (mo- 
menta per charge) of most of these particles (< 1 GV) do 
not allow them to escape into interplanetary space. 

In Fig. [TT] we show the emergent photon spectra gen- 
erated by re-entrant particles in comparison with the pre- 
viously computed spectrum formed during the primary 
proton-induced cascade. It can be seen that the spectral 
shape of the secondary components produced by re-entrant 
electrons and positrons in the 25-300 keV range is not sig- 
nificantly different from that of the primary spectrum and 
that re-entrant particles contribute together ~ 8% to the 
hard X-ray flux. Morever, we verified that this amount de- 
pends very little on the primary proton energy. We take this 
minor additional component into account when making our 
flnal predictions in ^ 

We also explored the contribution of re-entrant albedo 
particles generated by primary cosmic electrons. The corre- 
sponding additonal hard X-ray flux proves to be less than 1% 
of the flux generated by the primary electrons, i.e. < 10~^ of 
the flux due to cosmic protons, and thus can be neglected. 



3 PREDICTIONS FOR OBSERVATIONS 

In the preceeding section we demonstrated that cosmic par- 
ticles of different types and energies generate in the Earth's 



atmosphere hard X-ray radiation with a nearly universal 
spectrum, which can be approximately described by equa- 
tion ([T]). We also determined the dependence of the emer- 
gent flux on the primary particle energy. On the other hand, 
the spectrum and chemical composition of Galactic cosmic 
rays impmgmg on the Eartllfl are known to a good accu- 
racy, in particular thanks to the many experiments of the 
past decade. The spectrum is stable at high energies (above 
~ 10 GeV per nucleon) but varies significantly at lower en- 
ergies over the 11-year solar cycle due to modulation by the 
interplanetary magnetic field and solar wind, as well as with 
position on the Earth due to a rigidity cutoff imposed by the 
terrestrial magnetosphere. Since these effects are fairly well 
understood, we can combine the available information about 
the cosmic-ray spectrum with the results of our simulations 
for monoenergetic particles (^ to make predictions for ob- 
servations of the Earth's hard X-ray radiation by satellite- 
born instruments. 

3.1 Integrating over the cosmic-ray spectrum 

To describe the spectrum of cosmic protons as well as its 
dependence on the solar modulation potential (Zeff) corre- 
sponds to the average energy loss of cosmic rays inside the 
heliosphere, where Z is the charge number of a particle and 
e is the electron charge) , we use the approximation provided 

Note that after solar flares the Earth is bombarded by addi- 
tional cosmic rays with variable spectrum. This transient compo- 
nent is not considered here. 



8 S. Sazonov et al. 



bv lUsoskin et alj (|2005l ) , which is based on a summary of re- 
cent measurements. Specifically, we adopt for the differential 
spectrum of protons at 1 AU 

^ 1.9[{E + mE + ^ + 1.876)]-'-'' 

' ' l + 0.4866[(£ + $)(£ + $ + 1.876)]-i-255 

£:(g+ 1.876) 
+ + $ + 1.876) 
particles (s cm^ sr GeV/nucleon)-\ (5) 

where E (GeV) is the kinetic energy per nucleon and 
$ (GeV) — {Ze/A)<j> (here A is the atomic mass), so that 
for protons E = Ep and $ = e<f>. We refer the reader to the 
above-mentioned paper (see also e.g. iGaisser et al.ll200lD for 
references to the original data from the experiments AMS, 
BESS, CAPRICE, IMAX, MASS, etc, most of which agree 
with each other within ~ 10% at energies below ~ 200 GeV. 

The Galactic cosmic-ray spectrum is additionally mod- 
ified by the geomagnetic field. To a first approximation, 
adopted below, this modulation can be represented by a 
sharp cutoff below a certain value iicut (GV) of the mag- 
netic rigidity of a charged particle, whose dependence on 
the geomagnetic coordinates is well-known. In H3.3I below 
we estimate the systematic error associated with this ap- 
proximation. 

Proceeding with our analysis, we integrated the ap- 
proximation for Cp given by equation ([2]) over the spec- 
trum of cosmic protons given by equation ^ for speci- 
fied values of the solar modulation potential <j) and geomag- 
netic cutoff rigidity iicut. The integration was done up to 
Ep = 1000 GeV, protons with energies higher than 250 GeV 
contributing ~ 3% to the emergent radiation flux for the 
highest expected values of iicut ~ 15 GV and less for lower 
cutoff values. 

In Fig. [12] we show the predicted hard X-ray flux in- 
duced by cosmic protons as a function of iicut for two values 
of (j>: 0.5 and 1.5 GV, which are typical for solar minima and 
maxima, respectively (when the suppression of the low en- 
ergy cosmic-ray flux is weakest and strongest, respectively). 
One can see that the emergent hard X-ray flux is almost 
independent of iicut when iicut <^ 1 GV, as could be ex- 
pected since the flux of incident protons with rigidities be- 
low ~ 1 GV is strongly suppressed by solar modulation. As 
follows from Fig. 1121 the emergent hard X-ray flux is ex- 
pected to increase by a factor of ~ 2 during solar minima 
compared to solar maxima in regions with low geomagnetic 
cutoff rigidity, i.e. near the magnetic poles of the Earth. The 
solar-cycle variations become progressively smaller with in- 
creasing cutoff rigidity and reduce to ~ 10% near the geo- 
magnetic equator, where iicut ^ 15 GV. 

We next performed a similar calculation for Q-particles. 
In this case the Galactic spectrum is very similar to that 
of protons, while He composes ~ 5% (in particle number) 
of cosmic rays. The solar-modulated spectrum of a-particles 
can t hus be well approxim ated by equation ([5]) multiplied by 
0.05 l|Usoskin et al.ll2005h . Note that since for a given parti- 
cle energy the rigidity of a-particles is a factor of ~ 2 smaller 
than that of protons, the effects of the solar and geomagnetic 
modulation are accordingly modifled. The results of ^ in- 
dicate that the radiation flux produced by an a-particle of 
given energy is nearly equal to that induced by a proton 
of the same energy (at least at energies below 40 GeV). 



This makes the corresponding calculation straightforward. 
In Fig. [12] we show as a function of iicut (geomagnetic cut- 
off rigidity for a proton) the predicted hard X-ray surface 
brightness of the Earth's atmosphere due to bombardment 
by a-particles during solar minimum and maximum. The 
hard X-ray flux induced by a-particles proves to be ~ 22% 
and 28% at iicut < 1 GV and iicut > 1 GV, respectively, 
of that due to protons. 

We can estimate in the same manner the contribution of 
heavier cosmic particles. To this end we assume that atoms 
of the groups CNO, Ne-S and Fe have the same Galactic 
spectra as a-particles and protons, while their differential 
fluxes (per GeV/nu cleon) are 6.8%, 2.6% and 0.54% of that 
of a-particles (e.g. iFiorentini et al.ll200ll ). Since for these 
groups of atoms the relation Z = A/2 approximately holds, 
the helio- and geomagnetic effects for them are similar to the 
case of He. Assuming that, similarly to the case of He, the 
hard X-ray flux induced by a heavy particle is the same as 
that produced by a proton of equal energy, the hard X-ray 
surface brightness of the Earth's atmosphere due to bom- 
bardment by heavy cosmic particles can be readily estimated 

We can next estimate the hard X-ray surface brightness 
of the atmosphere due to bombardment by cosmic electrons 
and positrons. In this case we should use equation ([4]) in- 
stead of equation ([2]|. Also the magnetic rigidity effects are 
somewhat different from the case of protons and nuclei be- 
cause the rest energy of electrons and positrons is negligible. 
Adopt ing for the solar-mo dulated spectrum of cosmic elec- 
trons (|Webber et al.lll983l ) 

Jc(£;c, (t>) = 0.07(£c+e</>)-='-^ f TT-^^l (s cm'' sr GeV)-\ 

\Ec + e(f) J 

we obtain the dependence of the emergent hard X-ray flux 
on iicut shown in Fig. 1121 

Finally, the contribution of re-entrant albedo electrons 
and positrons can be roughly estimated by multiplying the 
hard X-ray flux due to cosmic protons and nuclei by 0.08 
(see H2.4p . since albedo particles are expected to re-enter 
the atmosphere in the opposite geomagnetic hemisphere in 
regions with similar geomagnetic cutoff rigidity to their es- 
cape positions. This emission component as well as the total 
atmospheric signal are shown as functions of iicut in Fig. 1121 



3.2 Approximate formula for the hard X-ray 

surface brightness of the Earth's atmosphere 

We found that over the region of practical interest: 0.3 < 
(fiiGV) < 2 and iicut (GV) < 50 - the dependence of the at- 
mosphere's hard X-ray surface brightness on the solar mod- 
ulation potential (f>, geomagnetic cutoff rigidity iicut and 
zenith angle 6 — cosO) can be well described by the 
following formula: 

^ _ 3m(1 + m) 1-47 X 0.0178 [(.^/2.8)°-^ + (0/2.8)^'^] ~' 
{l-F[iicut/(1.3<;!.o-25(l-K2.50O-4))]2}°-' 

(s cm^ sr)"""". (7) 

To calculate the emergent hard X-ray spectrum (from 1 cm^ 
of the Earth's atmosphere into a unit solid angle) for given 
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Figure 12. Top: Predicted hard X-ray surface brightness (inte- 
grated over all zenith angles) of the Earth's atmosphere during 
solar minimum {(f) = 500 MV) as a function of the geomagnetic 
cutoff rigidity (for a proton) iJcut • The thick solid line shows the 
expected total signal - coefficient C in the spectral approxima- 
tion given by equation jl}, which consists of several components 
(thin solid lines): due to protons, a-particles, CNO nuclei, Ne-S 
nuclei, Fe nuclei, cosmic electrons and positrons, and re-entrant 
electrons and positrons. The cumulative flux is well fit by equa- 
tion lO - the dashed line nearly coincident with the thick solid 
line. Bottom: The same but for solar maximum (0 = 1.5 GV). 
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Figure 13. Simulated spectra (solid lines) of Earth atmospheric 
emission induced by cosmic protons for = 0.5 GV, iJcut = 
0.5 GV and ip = 1.5 GV, iJcut = 15 GV (upper and lower curves, 
respectively). The dashed lines show the corresponding predic- 
tions for the 25—300 keV energy range based on equations JlJ and 
(O (integrated over all zenith angles). The dotted lines are sim- 
ulations based on the incident direction-dependent cutoff rigidity 
approach in the dipole approximation for the geomagnetic field 
(Rcut = 0.5 GV and i?cut = 15 GV approximately correspond to 
Am = 64.6° and 0°, respectively). All the simulated spectra were 
multiplied by 1.47 to account for the non-proton cosmic rays and 
re-entrant albedo particles. 



ing the outgoing radiation over all zenith angles. It can be 
seen that the computed spectra are in very good agreement 
with our analytic prediction in the 25-300 keV energy range. 
This confirmes the self-consistency of our treatment and the 
stated few per cent statistical accuracy of our simulations. 



values of (f>, Rcut and fi, one just needs to substitute the 
amplitude C given by equation ([7]) into equation ([l}. 

We first derived the above dependence for protons only 
and then proved that the combined contribution of the other 
cosmic-ray constituents, despite their somewhat different de- 
pendences on Rcut (see Fig. I12p . can be well represented by 
the additional factor 1.47 appearing in equation ((TJ. This 
factor also takes into account the estimated ~ 8% contribu- 
tion of re-entrant albedo particles. 

To verify our prediction for the atmospheric hard X-ray 
spectrum given by equations ([1} and ((Tjl , we performed a se- 
ries of simulations where the atmosphere was bombarded by 
the cosmic proton spectrum expected for specified (</>, i?cut) 
values according to equation ([5}. In Fig. [13] we present two 
extreme examples of the resulting emission spectra - one 
corresponding to observing a geomagnetic pole during so- 
lar minimum and the other to observing the geomagnetic 
equator during solar maximum, in both cases on integrat- 



3.3 Effect of the cutoff rigidity's zenith angle 
dependence 

In the preceeding discusssion we assumed that at a given 
position on the Earth the geomagnetic cutoff rigidity is the 
same for all directions of incidence. In reality, the effects of 
the geomagnetic field can be better represented by a cutoff 
rigidity that also depends on the zenith (0) and azimuthal 
{(,) angles. If the geomagnetic field is approximated by a 
dipole, t hen the cutoff rigid ity is given by Stormer's formula 
(see e.g. iLipari et al.lll99^ ): 

o ^ COS"^ Am , . 

2r2 [l + (l + cos3AmSinesinO'/2]2' 

where ~ 8.1 x lO'^^ G cm'' is the magnetic dipole mo- 
ment, r is the distance from the dipole center and Am is the 
magnetic latitude. In the atmosphere, i.e. near the Earth 

surface (r « r®), MQ/{2r^) ~ 59.4 GV, so at the geomag- 
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netic equator (Am = 0) i?cut takes values between 10.2 and 
59.4 GV, the vertical {9 = 0) cutoff being 14.9 GV; at the 
poles i?cut = for all directions. 

We ran a series of simulations where for given values of 
the solar modulation potential (j) and geomagnetic latitude 
Am the incident cosmic-proton spectra were dependent on 
the zenith and azimuthal angles, namely they were cut off 
at iZcut according to equation We then compared the 
produced spectra of atmospheric hard X-ray emission with 
those simulated assuming that regardless of the incident di- 
rection the proton spectrum is truncated at the vertical cut- 
off value: i?cut(Am, 6', = i?cut(Am,^ = 0). Both kinds of 
spectra turned out to be mutually consistent within a few 
per cent (see Fig. I13p, even at the geomagnetic equator. We 
note however that Stormer's formula generally underesti- 
mates the cutoffs because it allows penetration of charged 
particles with trajecto ries that would have intersected the 
surface of the Earth (| Smart fc Shea|[200i ). The emergent 
radiation flux may thus be smaller than predicted here by 
a few per cent. More accurate simulations should be based 
on a particle backtracking approach and use more realistic 
models of the geomagnetic field. 




20 40 60 



3.4 Observations of the Earth from space by a 
wide-field instrument 

Equation ((Tjl allows one to predict the hard X-ray surface 
brightness of the atmosphere or its part when observed from 
a spacecraft, provided that its orbit is higher than ~ 50- 
100 km. As a practical example, we may use this formula 
to predict the hard X-ray flux that would be measured by 
a satellite-born instrument whose field of view covers the 
entire terrestrial disk, which was the case in the recent IN- 
TEGRAL observations. In the dipole approximation for the 
Earth's magnetosphere, the spectral coefficient C in equa- 
tion (Ul) can be calculated as a function of the spacecraft 
geomagnetic latitude Am and altitude D (Fig. I14p . One can 
see that in the case of a low orbit {D — 400 km in our ex- 
ample), during solar minimum the fiux detected during the 
passage of one of the magnetic poles (Am) will be ~ 5 times 
higher than that measured during the passage of the equa- 
tor. In the case of intermediate and high orbits (20,000 and 
200,000 km in our examples) the relative amplitude of vari- 
ations is expected to be ~ 2. 

We point out that since the geomagnetic field is not ex- 
actly dipolar, in analyzing the data from a real experiment 
one should integrate equation (O over the geomagnetic rigid- 
ity map corresponding to a given observation. 



4 HARD X-RAY RADIATION FROM THE 
MOON AND PLANETS 

Obviously, the same mechanism of hard X-ray production 
should work in the crust of the Moon, so it is interesting to 
consider this case as well. Since the Moon has a negligible 
magnetic field, in its case both the incident cosmic-ray flux 
and the emergent hard X-ray flux are expected to depend 
on the solar cycle phase only. 

We performed simulations assuming that the Moon's 
crust has a density of 3 g cm~^ and consists of O, Si, Al, Fe, 
Ca, Mg, and Ti with mass fractions of 43, 21, 10, 10, 9, 5 and 



Figure 14. Hard X-ray spectral flux - coefficient C in equa- 
tion JTJ — that would be detected from the whole Earth by a 
wide-field instrument in different circular orbits (see notation on 
the figure), as a function of the spacecraft geomagnetic latitude. 
The thick and thin lines correspond to observations during solar 
minimum (</> = 0.5 GV) and solar maximum (</) = 1.5 GV), re- 
spectively. Note that the dependencies were rescaled as indicated 
along the vertical axis. 



2%, respectively. It turns out (see Fig.ll5[) that the spectrum 
emitted above ~ 100 keV is similar to that produced in 
the Earth's atmosphere and only the cutoff below 100 keV 
is stronger due to photoabsorption on heavy atoms. The 
Moon's spectral surface brightness can be approximated in 
the energy range 40-300 keV as follows: 



dNj _ 3/i(l + fi) 



dE 



5tt 



X 0.013 



V88 keVj 



\2.8j \2.8j 



E 

(s cm'^ sr keV)~^. 



(9) 



In writing down this formula we assumed, by analogy with 
the Earth's atmosphere, that the total radiation flux pro- 
duced in the Moon's crust is that due to primary cosmic 
protons times 1.4. Note that the corresponding factor for 
the Earth's atmosphere (1.47) included the small contribu- 
tion of re-entrant albedo particles. 

Similar results are expected also for Mars and Mercury, 
since the chemical composition of their solid crusts is sim- 
ilar to that of the Moon and they too have thin gaseous 
atmospheres and weak magnetospheres. Also the effects of 
interplanetary magnetic flelds on the flux of cosmic rays are 
not expected to be signiflcantly different in the orbits of 
Earth, Mars and Mercury. 

Interestingly, the Moon should appear as a fairly bright 
hard X-ray source when observed from a satellite orbiting 
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Figure 15. Expected spectra of cosmic ray-induced emission 
from the Moon's crust (integrated over all zenith angles) during 
solar minimum (</> = 0.5 GV) and maximum (</> = 1.5 GV) - up- 
per and lower histogram, respectively. The dashed lines show the 
approximation in the 25-300 keV energy range given by equa- 
tion (O. For comparison shown are the corresponding spectra 
emitted near the magnetic poles of the Earth {-Rcut = 0). 



the Earth. Indeed, on integrating equation Q over the Moon 
disk one finds that a detector located at distance D from the 
Moon will measure a spectral flux 

2 



,{D) = (4^8)xl0" 



3 X 10^ km 
D 



E 



keV 



+ 



keV 



photon 
s cm^ keV ' 



(10) 



where the lower and upper values correspond to observations 
during solar maximum and minimum, respectively. There- 
fore, a flux of '-^ 1 mCrab is expected at 100-200 keV. 

One can expect hard X-ray production by cosmic rays 
to be much less efficient on giant planets such as Jupiter and 
Saturn, since their typical cutoff rigidities (roughly propor- 
tional to the product of the magnetic field strength on the 
surface of the planet and its radius) are much higher than 
in the terrestrial case. Similarly, the spectrum of Galactic 
cosmic rays impinging on the Sun is truncated at TeV en- 
ergies jSeckel et al.lll99lh , and this strongly suppresses the 
generation of hard X-rays and gamma-rays in the solar at- 
mosphere. 



5 DISCUSSION AND CONCLUSIONS 

We have revisited the problem of the X-ray and gamma-ray 
emissivity of the Earth's atmosphere taking advantage of the 



recent high-quality measurements of the Galactic cosmic-ray 
spectra and the availability of efficient Monte Carlo codes in- 
corporating state-of-the-art descriptions of hadronic physics 
and electromagentic cascades. We found that emission gen- 
erated in the energy band 25-300 keV has a practically in- 
variant spectral shape and a flux that is expected to show 
signiflcant variations both with time (over the solar cycle) 
and with position on the Earth (due to geomagnetic effects) . 
We provided two simple expressions - equations ([TJ and (0 
- which allow one to predict the spectral intensity of at- 
mospheric hard X-ray emission for a particular observation 
from a spacecraft. It should be noted that our results are not 
directly applicable to observations performed from balloons 
floating at altitudes below ~ 50 km, since then in addition 
to the upward flux there will also be a signiflcant amount 
of downward radiation generated in the several overlying 
g cm~^ of air. 

The statistical accuracy of our predictions is better than 
5%, while the systematic uncertainty is probably dominated 
by the current uncertainty in our knowledge of the cosmic- 
ray flux, which is ~ 10%. We can try to assess the over- 
all accuracy of our model through direct comparison with 
real observations of the Earth from spacecraft. Among the 
few published measurements of atmospheric emission in the 
hard X-ray domain are those done in the 1970 's from the 
satellites Kosmos 461 jGolenetskii et al.l IWT^ ) and 1972- 
076B ((Imhof et al.iil97Q V The spectra obtained in the energy 
range ~ 50-1000 keV could be well described by a power law 
with a photon index {dNj/dE oc E~'^) T = 1.55±0.10 in the 
former experiment and F ~ 1.31-1.47 in the latter. These 
values are in good agreement with our calculated value of 
1.4 for the spect r al slo pe above ~ 50 keV. Furthe rmore, 
iGolenetskii et all (|l974l ') (see also lMazets et al.|[l975l ) mea- 
sured the dependence of the hard X-ray flux on geomagnetic 
cutoff rigidity over the range 3 < -Rcut (GV) < 17.5, and that 
result too appears to be consistent with our computations. 

Recently during a scanning observation of the Earth, 
the INTEGRAL observatory measured the spectrum of 
the cosmic X-ray background together with that of at- 
mospheric hard X-ray emission. Our analysis of the data 
|Ch urazov et al.l l2007f) showed that the atmospheric spec- 
trum at energies from ~ 40 to ~ 200 keV is consistent 
with the prediction of the present study, although the spec- 
tral shape was not well constrained by the data. When we 
adopted the spectral shape given by equation ((TJ, we found 
the best-flt amplitude of the atmospheric component to lie 
within 10% of (actually only ~ 4% higher than) the flux 
predicted by equation (0 for the relevant value of the solar 
modulation potential {(p — 0.45 GV) and taking into ac- 
count the geomagnetic cutoff r idigity map correspon ding to 
the INTEGRAL observations (|Churazov et al.ll2007l '). 

This good agreement of our model with observations 
suggests that the Earth could be used as an absolute cali- 
brator of future satellite-born sensitive hard X-ray detectors, 
complementing the Crab nebula. This possibility as well as 
the predictive power of our model could be tested soon via 
new INTEGRAL observations of the Earth performed at a 
different phase of the solar cycle. 

We also demonstrated that under the impact of cos- 
mic rays the solid surfaces of the Moon, Mars and Mer- 
cury should be about as bright in hard X-rays as (in fact 
somewhat brighter than) the Earth's atmosphere, making 
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Figure 16. Top; Total X-ray spectral surface brightness of the 
Earth (soild lines) — sum of cosmic ray-induced atmospheric emis- 
sion (dotted lines) and the reflected CXB radiation (dashed lines) . 
The two sets of curves correspond to observations of a geomag- 
netic polar region (Rcut = 0) during solar minimum (</< = 0.5 GV) 
and of the geomagnetic equator (iJcut = 15 GV) during solar 
maximum (0 = 1.5 GV). Bottom: The same but for the Moon. 
The two sets of curves correspond to observations performed dur- 
ing solar minimum and solar maximum. 

these objects potentially interesting for observations. It is 
important to keep in mind though that in a real experiment 
reflected CXB radiation will significantly contribute to the 
observed signal at energies below ^ 5 keV. Based on the 
computations of lChurazov et all |20o3), we show in Fig. [T6l 
several examples of hard X-ray spectra that could be ob- 
served from the Earth and Moon. This figure indicates that 
the cosmic ray-induced hard X-ray continuum could comple- 
ment the CXB-induced fluorescent X-ray lines in probing the 
chemical composition of the crusts of planets and the Moon, 
i.e. during the Bepi Colombo mission to MercurjQ. 
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